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ABSTRACT 

The opening mode stress Intensity factor and the associated crack 
mouth displacement are comprehensively treated using planar boundary 
collocation results supplemented by end point values from the litera- 
ture. Data are expressed in terms of dimensionless coefficients of con- 
venient form which are each functions of two dimensionless parameters, 
the relative crack length, and a load combination parameter which 
uniquely characterizes all possible combinations of tension or compres- 
sion with bending or counterbending Accurate interpolation expressions 
are provided which cover the entire ranges of both parameters. Applica- 
tion is limited to specimens with ratios of effective ha.f-heighc to 
width not less than unity. 


LIST OF SYMBOLS 

a crack length (or depth); see Fig. 1 

B specimen thickness; see Fig. 1 

c^.Cj numerical coefficients in E<r>. (11) 

E Younj's modulus 

E' effective modulus; see Discussion 

H effective half-height of specimen 

K stress intensity factor (Mode I in this context) 

M moment of complementary couple 

P resultant applied force 

s distance of knife edge from crack mouth (see Crack Mouth Dis- 

placement Slopes) 

v crack mouth displacement 

W specimen width; see Fig. 1 
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a relative crack length, a/W 

dimensionless stress intensity coefficient 
A dimensionless crack mouth displacement coefficient 

displacement combination parameter, arctanfv^/Vp) 
v Poisson's ratio 

o stress 

0 load combination parameter, arctan ( w 

Subscripts: 

1 opening mode of crack tip deformation 

M value for net section bending (u • 1) 

P value for net section tension (w * 0) 


1. INTRODUCTION 

This paper is a comprehensive treatment of the opening mode stress in- 
tensity factor, Kj, and the associated crack mouth displacement, v, for 
quasiplanar rectangular specimens with single side cracks of any uniform 
depth, loaded in any manner that can be represented sufficiently well by 
a linear distribution of normal forces across each end boundary. To 
justify such a representation of simple practical methods of load appli- 
cation, the ratio of effective half-height to width of the specimen 
should not be less than unity, preferably substantially greater for 
higher accuracy of application of the displacement data. The treatment 
of shorter specimens would require more refined, nonlinear modeling of 
the actual loading force distribution so that the interaction with the 
crack stress field would be adequately approximated. 

Results for certain cases covered here have been published previ- 
ously by the present authors, such as in Ref. [1], and by others cited 
in Rei . [2]. However, these previous results are not sufficient for a 
comprehensive treatment of the subject. The results of this study con- 
stitute a homogeneous body of data which were subjected to linear re- 
gression analysis to estimate consistency and precision. 

In this treatment the distribution of normal forces across an end 
boundary is characterized by the statically equivalent combination of 
resultant force, chosen to act through mid net section, and the comple- 
mentary couple. Either the force, or the couple, or both, can act in a 
positive sense so as to cause the crack to open in the vicinity of the 
tip (tension or bending), or can act in a negative sense to cause the 
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crack to tend to close near the tip (compression or counterbending). 

For some load combinations the effect on the crack mouth is the same as 
at the tip; for others, such as combined tension and counterbending, it 
may be opposite. In general, the force and the couple are to be regarded 
as Independent variables. Commonly, however, in practical testing, when 
force is applied by a single actuator, there is a proportional relation 
between the force and the couple which depends on the line of action of 
the force and on the relative crack length. 

The data are given in the form of dimensionless coefficients, r(K^) 
and A(v), which ate both functions of two dimensionless parameters, the 
relative crack length, a, and the load combination parameter, u>, equal to 
arctan(6M/P(W - a)), where P is the magnitude of the resultant force, 

M that of the moment of the complementary couple, U is the specimen 
width, and a is the crack length. Both T and A are linear functions 
of the quantity tan(n/4 - ui) , and the coefficients of the linear equa- 
tions aie finitely bounded functions of a which are expressed as 
rational algebraic functions for interpolation purposes. The ratio H/U, 
of the half -height to width of the specimen is the third Independent 
parameter which would be expected to have an appreciable effect on the 
values of r and A if it were less than unity. For the values con- 
sidered: 1, 2, and 4, however, the effect on T was found to be prac- 
tically negligible, and the eflect c* A erratic and of the order of 
5 percent at the most. 

The primary data for values of the relative crack length a other 
than the end points, 0 and 1, were obtained by the boundary collocation 
method of analysis as described in Ref. (3); the end point values were 
obtained from Refs. [2] or [4] except in one case which is discussed 
later. 


The purpose of this study is to serve as a prerequisite for the 
development of new test methods for measurement of such properties of 
materials as fracture toughness, fatigue crack propagation resistance 
and stress corrosion crack propagation resistance. The potential advan- 
tages of the full range of load combinations has by no means been thor- 
oughly explored. Nor has the possibility that the property oi Interest 
might be significantly dependent on the mode of loading. This possi- 
bility follows from the fact that the crack tip plastic zone character- 
istics will certainly be affected. 

In the application of results of planar analysis it is useful to 
distinguish between two types of quasi-planar practical plate specimens, 
as shown in cress section in Fig 1: (a) edge-cracked, and (b) iace- 

cracked In either case the crack length, or depth, must be practically 
uniform if the planar analysis is to apply ac-urately. Both kinds of 
specimen are necessary so that teats can be conducted with cracks ori- 
ented in different directions with respect to material texture; in par- 
ticular, in the direction that cracks have been found to occur in service 
of a structural member. One particular use of the face-cracked specimen 
follows from regarding it as the extr> me cast of the widely used part- 
through cracked, or surface cracked t pe of specimen. Fig. 1(c). The 
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analysis ot part -through cracked specimens Is essentially three- 
dimensional, and presents difficulties whl<h are not yet satisfactorily 
resolved [5], The present results provtde a bound on the three- 
dimensional analysis, and are therefore a useful guide as to Its accu- 
racy. Moreover, tests of face-cracked ipeclmem can provide a useful 
check on the reliability ot tests ot part-through cracked specimens, or 
serve to substitute tor them when convenient 


2. BOUNDARY C0L10CAT10N ANALYSIS 

Boundary collocation analyses were conducted for two dlflerent, but 
equivalent, types oi end boundary conditions, as shown in Fig. 2(a). The 
first tvpe Is a linear dl6trlbuticn of normal tractions across the end 
boundary, as shown at the top of the figure. The second type consists of 
unfform df strlbut lens of shear stress along pt olongat ions of the side 
boundaries, as shown at the bottom ot the figure. Both distributions are 
statically equivalent to the combination ot resultant iotwe, P, which 
acts througn mid net section, and complementary couple ot moment M. The 
use of these two distinct types ot boundary conditions provided a check 
ot the consistency of the method. It was found that the same results 
could be obtained with either type, but that the analysis had to be car- 
ried considerably further to reach satisfactory convergence when the sec- 
ond (shear) type of boundary condition was used. Satisfactory convergence 
is approached by successive runs with increasing numbers of boundary sta- 
tions until there Is no slgnlfican' variation In the result. The complete 
set of results presented here was therefore obtained with the first (nor- 
mal streas) type ot boundary conditions, and a substantial number of 
these results were confirmed with the second type. The purpose ot 
Fig. 2(b) is to show other equivalent expressions of the load combination: 
a linear distribution of normal tractions across the net section; and a 
pair of parallel forces acting at the corners. The**? expressions serve 
conceptual purposes; they are cf no concern to the boundary collocation 
analysis . 

The boundary collocation procedure used has been described In 
Ref. f 3 ] , and in gteater detail In Ref. 1 6 J . Briefly, a stress function 
which satisfies the ne' essary conditions on the crack boundary surface 
takes the form of an infinite series with unknown coef 1 icients . This 
series is truncated and the remaining coettlclents are determined from 
the conditions *hat the stress function and its derivative must satisfy 
the imposed bounds conditions at a finite number ol selected boundary 
stations. The derivation ot the boundary values of the stress function 
and its derivative lor the present cases are given in Appendix A. The 
number of undetermined coefficients can be equal to or less than twice 
the number of boundary stations 11 it Is less, then an overdetermined 
system of simultaneous equations results, and this set of equations is 
solved by a least squares best fit technique. This procedure has the ad- 
vantage of introducing smoothing into the numerical analysis, and has 
been adopted as a normal routine The value of the stress intensity fac- 
tor is directly proportional tc the coefficient of the first teim of the 
series, which is singular and dominates the stress and displacement 
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fields in the vicinity of c lie crack tip. Satisfactory convergence of 
this coefficient is sufficient to obtain the stress intensity factor 
accurately. Displacements at points not close to ihe crack tip Involve 
all the coefficients, not equally but with successively diminishing im- 
portance. Consequently the analysis must le carried considerably fur- 
ther to obtain satisfactory values of displacements than to obtain the 
stress intensity factor. 


3. FORM AND USE OF RESULTS 

The values obtained by boundary collo cation o f the stress Intensity co- 
efficient, T • r( a ,u>) • K|/(Cp ♦ Oty)»a(l - a), are listed in Table 1; 
those of the crack mouth displacement coefficient, A * A(a,u) 

■ E'v/(op + c M )a, are similarly listed in Table 2. As mentioned in the 
Introduction, these forms of dimensionless coefficient were chosen be- 
cause of their particular suitability for interpolation over the entire 
ranges of the two principal dimensionless variables: a “ a/W, and 

ui ■ arctanfo^/cp) • arctan(6M/P( W - a)). The form of T is essentially 
due to Kolter [4] and |7], and that of A is the natural complement. 
Discussion of the properties and use of these coefficients now follows. 

The value of the stress intensity factor, Kj , for any combination 
of forces is obtained, like a component of stress, by direct superposi- 
tion. In the present case the value of Kj for the combination of re- 
sultant force, P, and couple of moment M is: 

K " K P * Si " <r p C p + r M °M )/a(l ' a) (1) 

where the subscript I has now been dropped on the urn. standing that 
the present context refers exclusively to the first or owning mode of 
crack extension, and where 

Op • P/B(W - a) is the component of normal net stress due to P 

Oj^ ■ 6M/B1W - a)^ is the component of normal net stress due to M 

Ip ■ Kp/Op»a ( 1 - a) is the stress Intensity coefficient for uniform 
normal net stress, as gi en in Table 1 for u> - 0 

" W* (1 “ * 8 t * ie 8tre8S Intensity coefficient for pure 

bending net stress distribution, as given in Table 1 for 
ui ■ ir/2 

The combined stress intensity coefficient is defined as: 


r • + o M Wa(l - a) 


( 2 ) 
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■ ( Vf + 0 m‘m ,/(0 P ♦ °M> 

(3) 

• (Tp ♦ ten w)/(l ♦ tan w) 

(4) 

where tan w • c^/o p • 6M/P(W - a) 


Similarly, for the crack mouth displacement, v: 


v - Vp ♦ v M - (Op Ap 4 o M A M )a/E' 

(5) 


where 


A p • E'Vp/Opa 

4 m • E 'W 

2 

E ■ E, Young's modulus, for piano stress, or E/(l - v ) for plane 
strain, where v is Poisson's ratio 

The effective value of E' for a practical specimen will lie between 
these limits, and will depend on the ratios a/W and B/W. It is best 
established by direct experiment on actual specimens or scaled-up models. 

Also, the combined crack mouth displacement coefficient is defined 

as : 


A ■ E'v/(o p ♦ o M )a (6) 

• <°p Ap ♦ o M A M )/(o p + o M > (7) 

■ (Ap ♦ A^ tan ui)/(l + tan ui) (8) 

The advantage of the form of the parameter u> is that every possible 
state of load combination corresponds to one, and only one, value of u> 
in the range 0 <. u> <. 2" Figure 3 shows how the four kinds of mixed 
load combinations might be obtained with a single actuator. Also, eight 
special cases can be distinguished as follows: 



tan ui 

State of load combination 

0 

0 

Simple net tension 

i»/4 

1 

Balanced tension and bending 

n/2 

OD 

Simple net bending 

3w/4 

-1 

Balanced bending and compression 

ir 

0 

Simple net compression 

5ir/4 

1 

Balanced compression and counterbending 

3n/2 

— no 

Simple counterbending 

7n/4 

-l 

Balanced counterbending and tension 
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Equations (4) and (8) can be expressed in linear form: 

r • (r p ♦ r M )/2 ♦ (1 - tan u)(r p - r M )/2(l ♦ tan u.) 


♦ ! m >/2 + 

Up - r M )tan("/4 - w)/2 

(9) 

* V' J * 

(Ap - A M )tan("/4 - w)/2 

(10) 


These linear equations enable the results given in Tables 1 and 2, for 
three values of w, to be subjected to linear regression analysis, which 
provides an objective measure of their precision. The sample size of 
three is minimal, but Eqns. (9) and (10) are true functional relations, 
not statistical ones, and the three values of the Independent variable, 
tan("/4 - ui) • 1, 0, and >1, are exact; consequently any variation about 
the regression line is attributable solely to scatter of the boundary 
collocation values of the dependent variable, F or A. 

The measure of precision calculated was the ratio of the standard 
error of the estimate to the mean of the three values of the dependent 
variable, T or A. The standard error of the estimate is the positive 
square root of that part of the variance of the dependent variable not 
accounted for by its regression on the independent variable. The preci- 
sion of the values of F in Table 1 is better than 0.1 percent in all 
cases. In Table 2 the precision of the values of A is also better than 
0.1 percent when H/W is 4, but when H/W is 2 or 1 about half of the 
results are less precise, the worst being hardly better than 1 percent. 
These less precise results are for a • 0.1, 0.2, 0.3, or 0.8, not for 
intermediate values, which reflects the fact that the boundary colloca- 
tion method is subject to greater round-off errors when the relative 
crack length is either small or large. 

Since both T and A are linear in the parameter tan(*/4 - ui) . 
they are linearly related to one another. By elimination of the common 
parameter in Eqns. (9) ard (10): 

I - ■ C 1 » c 2 A (ID 


where : 

2 cj(a) • p + r ^ ~ (Ap + A^) ( p - r^)/(Ap “ A^) 

c 2 (a) - (r p - r M )/(A p - A m ) 

Alternatively, Eqn. (11) can be inverted to give a as a linear function 
of r. Either equation can be used In a similar manner to Eqns. (9) and 
1 1 H) .'or linear regression analysis. The correlation coefficients for 
Eqn. (11) were better than 0.99999 in all cases, and the precisions of 
the regression lines were better than 0.1 percent for H/W * 4. It is. 
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therefore, apparent that (ha boundary collocation values of I and A 
are hlgh.y consistent 

There Is a further linear relation for the ratio 
I/A • Ka^2/g V (1 a )l/2 f which Is useful when the crack mouth displace- 
ment is the controlled variable In fatigue crack generation or crack ex- 
tension resistance experiments. Let v M /v p • tan 6 • (A M /A p )tan u. , then 
It is easily shown that: 

2(T/A) • C p M p ♦ r M /fi M ) ♦ ( r p/A p - Ig/A^umf"/* - 0) (12) 


4. END POINT VALUES AND INTERPOLATION EXPRESS'ONS 

The combined coeiflclents ' and A are exact functions of the 
parameter w by definition - Eqns . (4) and (8) - but the particular co- 
efficients in these equations, namely: Ip • 1(1,0), I ^ * I(a,-/2), 

Ap ■ A(a,0) and • A(a,-/2), are unknown functions of the parameter 
a. It is therefore desirable to provide suitable interpolation expres- 
sions for these coefficients. These expressions should agree exactly 
with the end point values of the coefficients at o • 0 and 1, which may 
be assumed to be known very accurately. They should also be acceptably 
consistent with the end point slopes, di/da and dA/da, though not so 
stringently as to impose undue complexity on the forms of the interpola- 
tion functions. Finally they should agree with the Intermediate values 
of the coefficients, and should not deviate systematically from those 
values. The intermediate values used In the study of interpolation ex- 
pressions were those obtained by linear regression analysis ot the bound- 
ary collocation results for H/W • 4, which therefore Involve the results 
for all three values ot w In Tables 1 and 2. 

The various end point values are given in Table 3, each with a 
symbol in parentheses to identify its source. Ten of these values were 
obtained from either Kef. (2] or (4). Estimates of the other six were 
obtained by extrapolation of iunctions of the variables such that the 
slopes of the plots were nearly constant near one or the other of the 
end points. These estimates are identified by the symbol (T) for text. 
One of them is the coefficient i(l,0), and Its value of -3.26 is esti- 
mated to be accurate within about 1 percent. The other five are end 
point slopes, and may be somewhat less accurate, but are correspondingly 
less important for the present purpose It was established that the 
order of dA/da at i • 0 is -1 10^2 for both values of u> from the 

fact that plots of A against approach a common s.ope of -5.4 as 

a approaches zero. 

Interpolation expressions of general conic form which cover the full 
range of a from 0 to 1 are discussed by Bentham and Kolter in Ref. [4] 
for the case of and also for the case of uniform tension applied re- 

motely across the gross section, which corresponds in present terms to 
the combined load case where w ■ arctan(3a/(l - a)). In each case these 
authors determined the six disposable ccefflcients of the conic form from 
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the values of r at the end ,*;ints and those of Its first and second de- 
rivatives, without recourse to any intermediate values. Since the values 
of these coefficients are not given In Rei . (4), the present authors de- 
termined them from the data ot that reference, but found that the result- 
ing expressions were not in suitlclently good agreement with the boundary 
collocation data A modified approach was then tried in which only the 
end point values of T and its first derivative were employed, the re- 
maining two conditions required to determine the six coefficients being 
obtained from the boundary collocation data. This approach was no more 
successful than the first, so it was concluded that the general conic 
form was not the most appropriate for the present purpose. 


Other forms of Interpolation functions are given by Tada et al. in 
Kef. (2) for both and A^, and for the corresponding coefficients 

for the combined load case of uniform gross tension. These functions, 
however, do not agree sufficiently well with the present boundary collo- 
cation results, though they do agree with the less extensive results of 
Kefs. [1] and |8|. Consequently it was considered necessary to explore 
other forms of interpolation expressions. The most satisfactory inter- 
polation expressions devised were the rational algebraic functions which 
follow: 

T p - 1.9387 - 1.468 a - 4.76 a(l - a)/(l ♦ a) 2 (13) 


r M • 1.9887 - 1.326 a - (3.49 - 0.68 a + 1.35 a 2 )a(l - a)/(l + a) 2 

(14) 

A p - 5.84 - 9.1 a - 4.1 a 1/2 (l - a) 3/2 /(l ♦ 0.1435 a) (15) 

& M • 5.84 - 2.3 a - 3.42 o 2 + 2.52 a 3 - 4.1 a I/2 (l - a) 3/2 /(l + 0.1435 a) 

(16) 

In each of these expressions the leading coefficient is the value of the 
function for a ■ 0, and the remaining coefficients are related by the 
value for a » 1 , so that each expression has four or less independent 
numerical coefficients to fit eight boundary collocation values. More- 
over, Eqns . (13) and (14) are devised to have features in common, as are 
Eqns. (15) and (16), so that the general interpolation expressions which 
correspond to Eqns. (4) and (8) each have only five independent numerical 
coefficients to fit twenty four boundary collocation values (Tables 1 
and 2 for H/W ■ 4). The interpolation expressions therefore involve 
considerable smoothing ot the boundary collocation results. For all six 
expressions the algebraic mean of the deviations of the boundary colloca- 
tion values is less than 0.1 percent, and the deviations are unsystematic. 
For Eqns. (4), (13), and (14) the individual deviations are within plus 
or minus 0.5 percent, except for a ■ 0.1 where the equations overesti- 
mate the boundary collocation values by about 1 percent. It is quite 
possible that these boundary collocation values tor such short cracks are 
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In error to that extent. For Eqne (8), (15), and (16) the individual 
deviatlone are all within plue or minua 1 percent. From theee varloue 
considerations it ie concluded that the Intercalation expressions given 
here are sufficiently precise and reliable for practical purposes. 


5. CRACK MOUTH DISPLACEMENT SLOPES 

It is sometimes convenient to use separate knife edges, which are at- 
tached to the specimen by screws or adhesive, tor mounting the clip gage 
which is used to measure the crack mouth displacement. In such cases it 
is necessary to know the values of the displacement coefficient A g at 
gone small distance, s, away from the crack mouth. For this reason 
values of the crack mouth opening slope dA/dx were obtained and are 
listed in Table 4 Th* distance x lb measured away from the specimen 
edge along the extens*' .i of the crack line. Then, at a position x ■ s, 
the value of the displacement coefficient A g • A 0 ♦ s dA/dx, where A 0 
is the value of A at the crack mouth as given in Tables 2 and 3, or by 
Eqns (8), (15), and (16). In Table 4 the values obtained by boundary 
collocation, for a * 0.1 to 0.8, are somewhat less accurate than those 
of the displacement coetflclent Itself in Table 2. The reason Is that to 
obtain the derivative Involves taking the dltierence between the displace- 
ment at the crack mouth and that at some small distance along the crack. 
The relative error of this difference is likely to be considerably greater 
than that of either displacement. Tc obtain the derivatives with substan- 
tially greater accuracy would require that the boundary collocation anal- 
ysis be carried considerably further. This was not considered to be war- 
ranted since it should not be necessary for the supplementary term 
s dA/dx to exceed about one tenth ot the principal term A 0 . As a prac- 
tical rule the distance s should not exceed one-tenth of the crack 
length. 

The value in Table 4 for a • 0 was obtained by extrapolation of 
the inverse square roots of the boundary collocation values to their 
common *'alue. The values for a ■ 1 are the same as the values of A 0 
because the crack tip and the "hinge point" are then both located at the 
back surtace of the specimen. 


6. DISCUSSION 

In Keis. ( 1 1 and |8| the present authors published stress Intensity and 
crack mouth displacement coefficients, respectively, tor side cracked 
specimens loaded in uniform tension across the gross section. This 
boundary condition is equivalent to tensile loading through frictionless 
pins located at mid width, as shown by the inset in Fig. 4, and it in- 
volves the relation: 0m/np * tan w * 3a/(l - a). In this figure the 

straight lines correspond to Eqn. (9) with coefficients obtained from the 
interpolation expressions, Eqns. (13) and (lk). The points marked by 
circles on this chart are boundary collocation values of T tor uniform 
gross tension which are derived either from Kei . |1] or from recent sup- 



n 


pi ement ary results. The maximum devlat Ion of any of these points from 
the corresponding line Is less than 1 percent Thus the earlier results 
for this special case ate in very good agreement with the Interpolation 
expressions presented here. 

Figure A also illustrates the point that the state of loading for 
center pin loaded specimens varies from net tension (w • 0) when a • 0 
to net bending (w ■ * / 2 ) when a • 1. At » • 1/A the value of w Is 
"fit, and the net bending stress is exactly equal to the net tension 
stress; for greater values of o the bending stress Is predominant. Evi- 
dently, a face cracked specimen loaded In this manner does not represent 
the extreme case of a part-through cracked specimen, as has sometimes 
been mistakenly assumed For this purpose the specimen would have to be 
loaded In such a manner as to maintain uniform net tension for all values 
of a. Alternatively, It might be more practical to load the specimen so 
as to maintain uniform displacement across the width at some appropriate 
distance from the rra'k These two alternatives are not exactly equiva- 
lent, of course, and might produe significantly different test results. 
This is a matter for further Investigation. 

Figute 5 is the chart of the crack mouth opening displacement coef- 
ficient £ which corresponds to Fig A. The straight lines represent 
Eqn (10) and the coefficients are obtained from the Interpolation ex- 
pressions, Eqns. (15) and (16). The circles represent boundary colloca- 
tion results tor uniform gross tension from Kef. |8|. The deviations of 
these earlier results from the lines are less than 1/2 percent, except at 
a • 0.7 which Is 2 percent below the corresponding line It is notable 
that these values ot A for uniform gross tension do not differ greatly 
from those for net bending (tan (*/A - u) • -1), but are In strong con- 
trast to the value* for uniform net tension (tan("/A - u>) • 1). Whereas 
the minimum value of £ for grosB section tension Is 2. A, at a • 0.7, 
that for net section tension is -1.26, at a » 1. 

The reason for the rapid decrease ot & with increase of a under 
net section tension Is that the equivalent combination of torce and 
couple causes c ount e rbending of the two parts of the specimen on either 
side of the crack, and this counterbendlng is greater the longer the 
crack In consequence, although the crack tip is always opened by this 
mode of loading (the stress Intensity is always positive), the crack 
mouth displacement is negative when a exceeds about 0.65 In practice 
the crack Is usually simulated by a pointed 6lot or note) from which ex- 
tends a short fatigue crack. The slot mouth displacement can therefore 
be negative It is important to appreciate, however, that the present 
results do not apply if there is contact along any part of the crack 
faces, because this would Induce compressive contact stresses which were 
not Imposed In the present boundary collocation analysis. The effect of 
such contact stresses can be analyzed by various methods if their magni- 
tude and distribution can be determined 

The ctack mouth displacement data are needed for several different 
purposes, such as in the procedure of ASTM Standard Method ot Test E 399 
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for determination of the Kj c plane strain fracture toughness, or 'it 
displacement control of fatigue cracking in an electrohydraullc closed 
loop loading system. Another Important purpose, less commonly appreci- 
ated, is for analysis of test results when the measured force applied to 
a cracked specimen Is complemented by an unmeasured bending moment. This 
can occur unintentionally if loading fixtures are designed with insuffi- 
cient care or appreciation of the nature of the stress distribution im- 
posed on the specimen by the f.xture It will occur when force is trans- 
mitted to a face cracked specimen by wedge or hydraulic grips In order to 
obtain a uniform d'strlbutlon of end displacement across the specimen 
width. 


The compliance of the specimen, vB/P, can be measured by recording 
P against v while the specimen is lightly loaded and unloaded. Then, 
providing that the crack length is accurately known and practically uni- 
form, the value of tan w and thus can be determined as follows; 

from Eqn. (6) : 

E'v . o p yq) ♦ t M (a ) t , y 

(o p ♦ o M )a ‘ " (Op ♦ o M ) (1 + tan u)o p a 


Hence : 


tan id 



(E'v/o p a) - A p 


(E'Bv/P)(l - o)/o - A p 


“p * w 

1 + tan id 


(17) 


(18) 


Since the value of a is known, those of Ap and A M can be ob- 
tained from Eqns . (15) and (16), respectively, vB/F is determined exper- 
imentally, so that it only remains to determine the value of E' in 
order to obtain the value of tan u> from Eqn. (18), and hence the value 
of o^. The quantity E' Is equal to E/(l - x^) , where E is Young's 
modulus and x is equal to 0 for a homogeneous state of generalized 
plane stress, and is equal to Poisson's ratio for a homogeneous state of 
plane strain. In practice the valu<: of x for most metallic naterials 
is between 0 and 0.3, depending upc n the three parameters: a, id, and 

B/(W - a), but independent of size scale. For a given specimen form and 
material the value of x can be determined from a series of compliance 
experiments in which the values of a and w are accurately known. It 
seems likely that a systematic series of such experiments with several 
different materials might reveal generalities about the dependence of x 
on the three dimensionless parameters, and on Poisson's ratio. The 
authors are not aware of the existence of such information at the present 
time. 
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?. SUMMARY 

Comprehensive lntcrm.it Ion is presented on the opening mode stresr, in- 
tensity factor, Kj, and on the associated crack mouth displacement, v, 
for quasl-planar rectangular specimens with single transverse side cracks 
of any uniform depth, loaded in any manner that can be represented suffi- 
ciently well by a linear distribution ot normal forces across each end 
boundary The data tor values of the relative crack length, a, other 
than the end points, 0 and 1, were obtained by boundary collocation 
analysis; the end point values were available from the literature, except 
in one case which was obtained by a reliable method of extrapolation. 

The information is given in the form of dimensionless coefficients, 

1 (Kj) and A(v), which are both tun.tlons of two principal dimensionless 
parameters, a, and the load combination parameter, w. Both I and & 
are linear functions ot the quantity tan( r /4 - ui) , and the coefficients 
of the lined' equations are finitely bounded functions of a. Rational 
algebraic interpolation expressions are given for these functions of a 
which fit the primary data within the precision of that data, and are 
considered to be accurate to within 1 percent or better. The effect on 
T and A ot the third dimensionless parameter, the ratio H/W ot half- 
height to width of rne specimen, is 6hown to be very weak providing that 
it exceeds unity. 

The purpose ot the information l c to serve ns a prerequisite for the 
development of new test methods for measurement of such properties of 
materials as fracture toughness, fatigue crack propagation resistance, 
and stress corrosion crark propagation resistance. The potential advan- 
tages of the full range ot load combinations (tension or compression com- 
bined with bending or counterbending) has by no means been thoroughly ex- 
plored. Nor has the possibility that the property of Interest mi'**''; be 
signif irant ly dependent on the mode cf loading This possibilit follows 
from the fact that the characteristics ot the crack tip plastic zone will 
certainly be affected 
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APPENDIX: STRESS FUNCTION AND STRESS FUNCTION BOUNDARY CONDITIONS FOR 

SIDE CRACKED PLATES 

The form of the stress function in terms of the polar coordinates, r and 
$ as shown in Fig. A, is 


><r 


n+1/2. 

L—J A 2i»-1 


cos(n - 3/2H 


- + ?)♦] 


n-1 


n+1 » r / 
r A 2n (cos(n - 


1)$ - cos(n + 1)$] 


(Al) 


Its application is discussed in detail in Ref. (3). The values of the 
stress function, x» and its normal derivatives, 3x/3x and 3x/3y , along 
the prescribed model boundaries (AB, BC, and CD in Fig. A.l and AB , BB' , 
B'C' , C'C, and CD in Fig. A. 2) are derived from known boundary tractions 
as shown in Fig. A, and are as follows: 

From the loa* condition of Fig. A.l where P ■ 0 o BW and M ■ BW^o oM /6. 
Along AB 


X • 0 


(A 2) 


& 

3x 


0 


(A3) 


Along BC 



+ 



2* - 
ay 


o 


Along CD 


X 


PW , M 
Tif + B 


(A4) 

(A5) 


(A6) 


2 * - l 

3x B 


(A7) 
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From the load conditions of Fig A. 2 where t c and 
formly distributed side shear tractions, we have the 


T u are the uni- 
followlng: 


Along Afi 

X • 0 



0 


Along BB' 

X • 0 



r c fy - H) 


(AH) 

(A9) 


(A10) 

(All) 


Along B‘C‘ 

X • t £ (V - H)(x + a) 


( A12) 



Along C'C 


X - t c (V - H) W 


(A13) 


(A14) 


4* - x (V - y) ♦ t (V - H) 
u ' c 

Along CD 

X - t c (/ - H) W 

• (r + T )(V - H) 

9x u c 


(A15) 


(A16) 


(A1 7) 
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Table 1 

Boundary collocation valuer of the stress Intensity coefficient T 
obtained with normal traction boundary conditions 


H/U 

U) 



For a 

(• a/W) values of: 



0.1 

0.2 

0. 3 

0.4 

0.5 

0.6 

0.7 

0.8 



r - k, 

/(o p ♦ 

^l/aTT 

- a) 



4 

0 

1.470 

1.164 

0.961 

0.822 

0.725 

0.659 

0.614 

0.581 


n/4 

1.528 

1.252 

1.064 

.930 

.832 

.759 

.704 

.662 


*/2 

1.585 

1.339 

1.167 

1.039 

.938 

.858 

.794 

.743 

2 

0 

1.466 

1.164 

0.961 

0.822 

0.725 

0.659 

0.613 

0.580 


■n/4 

1523 

1.251 

1.064 

.930 

.832 

.758 

.704 

.661 


n/2 

1.583 

1.339 

1.167 

1.038 

.938 

.858 

.794 

.742 

1 

0 

1.468 

1.163 

0.961 

0.822 

0.725 

0.659 

0.614 

0.580 


n /4 

1 ,526 

1.250 

1.065 

.932 

.832 

.758 

.704 

.661 


1/2 

1.582 

1.338 

1 168 

1.039 

.938 

.858 

.794 

.741 


Table 2 

Boundary collocation values of the crack mouth opening coefficient 
A obtained with normal traction boundary conditions 


H/W 

U) 



For a 

(» a/W) values of : 


0.1 

0.2 

03 

0.4 

0.5 

0.6 

0.7 0.8 



A 

• E'v/(o p ♦ o M 

)a 


4 

0 

3.800 

2.770 

1.855 

1.055 

0.331 

-0.358 

-1.035 -1.730 


n/4 

4.125 

3.405 

2.763 

2.223 

1.757 

1.335 

.934 .540 


n/2 

4.450 

4.040 

3.670 

3.390 

3.183 

3.027 

2.903 2.810 

2 

0 

3 690 

2.735 

1.850 

1 055 

0.331 

-0.357 

-1.037 -1.733 


n/4 

3.985 

3.321 

2.741 

2.220 

1 757 

1.335 

.933 .532 


n/2 

4.400 

3.960 

3.675 

3.381 

3.183 

3.027 

2.901 2.800 

1 

0 

3.860 

2.780 

1.855 

1.055 

0.327 

-0.367 

0 

H 

1 

o 

1 


n/4 

4.200 

3.420 

2.763 

2.225 

1.763 

1.337 

.929 .536 


n/2 

4.550 

4.060 

3.670 

3.394 

3.190 

3.032 

2.903 2.795 



18 


Tabic 3 

End point values of the coefficients and their derivatives 
Sources in parentheses: Kefs. [2] or (4] or text (T) 


U) 

a 

r 


df/do 

6 



dA/ da 


0 

0 

1.9887 

(4) 

-6.96 

(4) 

5.84 

(2) 

ORDER 

<-l/0 1/2 ) 

(T) 

0 

1 

.5204 

(2) 

-0.4 

(T) 

-3.26 

(T) 

-91 

(-l/Ql/2) 

(T) 

*/2 

0 

1.9887 

(4) 

-5.40 

(4) 

5.84 

(2) 

ORDER 

(T) 

i/2 

1 

.6629 

(4) 

-0.332 

(4) 

2.64 

(2) 

-1.6 


(T) 


Table 4 


Values of 
dA/d (x/W) 

the crack 

mouth opening 

slope 

a/W 

For 

u values of: 



0 

n/4 

n/2 


d£/d (x/W) - (E'W/oa)dv/dx 

0 

48 

48 

48 

0.1 

15.3 

18 2 

21.1 

.2 

6.6 

9.3 

12.0 

.3 

2.6 

5.3 

8.0 

.4 

.44 

3.2 

6.0 

.5 

-0.b4 

2.0 

4.9 

.6 

-1.7 

1.2 

4.2 

.7 

-2.3 

.63 

3.6 

.8 

-2.9 

.27 

3.4 


.9 
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l.gure 1. ( foss Mtl'um through . raU planes at m tanguiar 

plate vecmens »ith 141 edge cratk <bi late crack ic> part 
through semieiiiplical 1 raik Mol* that the isidth dimension 
IV IS always taken to be in the tame direction as the i-rack 
length 4. dnrt the truck nest H taken to tie nor mal to that 
direction. 
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figure ?. - Specimen model dimensions. equivalent end boundery condi 
lions end stress distributions, ui Boundery collocation model shotting 
normel tractions case at top arid shear tractions case at bottom. |>i 
Equivalent linear net stress distribution model shotting equivalent 
forte pair atbo*tom. 





O Boundary collocation values lor center pm 
loaded specimen, ret. [1] 

A [no point values tor center pin leaded specimen 


Bending tan m'4 -w> -lop o M i«Op »Oyi tension 

figure*. Stress intensity coefficient, r. versus the (unction 
tame i yi o ( the load combination paramder. » The straight 
lines are the results ot the present work tor u • 0. vii , and 
fill the superimposed points are from earlier work on the case 
tor which tan u • Tali I 01 . 







Bendrnq tam*4 J-tflp O^i'Op* d v ' tension 


hqureS. Crack mould displacet'ient coetticient A, versus 
tan i " 4 wl. The straight linn are the r esults at the present 
iw* the superimposed points are Irorr. ear net work on the case 
tor which tan u ■ to i| or 


NOR MAI TRACTIONS SHIAR TRACTIONS 



Figure A Alternative boundary collocation models 



